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We discuss quantum algorithms, based on the Bernstein- Vazirani algorithm, for finding which vari- 
ables a Boolean function depends on. There are 2" possible linear Boolean functions of n variables; 
given a linear Boolean function, the Bernstein- Vazirani quantum algorithm can deterministically 
identify which one of these Boolean functions we are given using just one single function query. 
The same quantum algorithm can also be used to learn which input variables other types of Boolean 
functions depend on, with a success probability that depends on the form of the Boolean function 
that is tested, but does not depend on the total number of input variables. We also outline a procedure 
to futher amplify the success probability, based on another quantum algorithm, the Grover search. 



1 Introduction 

In the oracle identification problem, we are given an oracle from a set of possible Boolean oracles, and 
our task is to determine which one we have CD-El. The complexity of the problem is measured by the 
number of times we must query the oracle in order to identify it. Both the Bernstein- Vazirani |4] [5) 
and Grover quantum algorithms (6) [71 solve this type of problem. The Bernstein- Vazirani algorithm 
identifies linear Boolean functions with a single function query, and Graver's search algorithm finds 
marked elements in a database with N elements using ff(y/N) queries. 

Consider the following task. We are given a black box that evaluates a Boolean function f{x\ , xi , . . . x n ) 
that maps {0, 1}" to {0, 1}. The function depends on the values of at most m of the variables and is in- 
dependent of the other n — m. Such a Boolean function is called a junta, and, if it depends on only 
one of the variables, it is called a dictatorship. Our task is to find which of the variables the function 
depends on. We shall show how a variant of the Bernstein- Vazirani algorithm can solve this problem. 
Recently, Rotteler presented a quantum algorithm for identifying quadratic Boolean functions [8]. Atici 
and Serviedo discuss a quantum algorithm for identifying ^-juntas, essentially based on the Bernstein- 
Vazirani oracle @. The quantum algorithm we outline is simpler; moreover, we also present a method 
to further increase the success probability, based on Graver's quantum search algorithm. 

The paper is arranged in the following way. In section [U we review the Bernstein- Vazirani algorithm. 
In Section |3l we show that this quantum algorithm can also be used for the more general task of finding 
variables other types of Boolean functions depend on. In section [4] we show how a method based on the 
Grover search can be used to further increase the success probability of finding variables the Boolean 
function depends on. We finish with Conclusions. 
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2 The Bernstein- Vazirani algorithm 

The Bernstein- Vazirani algorithm is a one-shot quantum algorithm HIH. It solves the following prob- 
lem. One has a black box that evaluates a linear Boolean function, given by 

n 

f(x) =y-x=Y,yjXj, (1) 

7=1 

where the addition is modulo 2 and y is a fixed, but unknown, n-bit string. We want to find y. The 
Bernstein- Vazirani algorithm does this with one evaluation of the function. It does so by mapping the 
functions to vectors in an Af-dimensional Hilbert space, M 1 = (8>"^2, where N = 2" and H2 is a two- 
dimensional Hilbert space. The computational basis vectors of J4?2 are |0) and |1), and the basis vectors 
of ffl are labeled by ra-bit strings \x) = \x\ ) tg) \%2 ) ■ ■ ■ <8> \x„ ) . The function y ■ x is mapped to the vector v y , 
where 

( x \v y ) = -L(-iy*. (2) 



N 

These vectors are orthonormal, i.e. (v y |vy) = S y y , and they constitute an orfhonormal basis of J4? known 
as the parity basis [4]. This follows from the identity 

£ ( D >v <>,„• (3) 

xe{0A}" 

Because the vectors are orthonormal, they are perfectly distinguishable, and so with one measurement 
we can perfectly determine which function the black box is evaluating. 

This is actually accomplished by using a circuit consisting of Hadamard gates and an /-controlled- 
NOT gate. The Hadamard gate is the unitary transform 

|0> ^(|0> + |1» 

|1> ^(10) - 11>). (4) 
If we apply n Hadamard gates, one to each qubit in the state \x), we obtain 

H**\x) = -L £ (-ir%), (5) 
V^V z e{o,i}" 

where, as before, we have set N = 2". The /-controlled-NOT gate, where / is a Boolean function, acts 
on n + 1 qubits in the following way 

U f \x)\z) = \x)\z + f(x)), (6) 

where \x) is an «-qubit computational basis sate, \z) is a one qubit state (z, = 0, 1), and the addition is 
modulo 2. Now, the input state to the Bernstein- Vazirani circuit is the (n + l)-qubit state 

|^v,) = -^|00...0)(|0)-|l)). (7) 
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We first apply n Hadamard gates, one to each of the first n qubits, and then the /-controlled-NOT gate, 
giving us 

|^> ^ -L £ (_1)/M| X )(|0)-|1)) (8) 



xe {Q.l}n 

Next, we again apply n Hadamard gates to the first n qubits yielding 

= ^TTj I I (-i) /w+M k>(|o>-|i». (9) 

"V z xe{o,i}" z e{o,i}" 

Discarding the last qubit (it is not entangled with the others, so this has no effect) and expressing this 
result in terms of the vectors |v y ), we find the «-qubit output state 

\Yout)= £ (v z |v/)|z), (10) 
ze{o,i}" 

where we have defined the vector vy to have the components (x\vf) = (l/y/N)(— l)^ x \ Now, if we know 
that f(x) is of the form f(x) = y-x, then we just get the vector \y) as our output, and when we measure 
I Wout ) in the computational basis, we find the «-bit string y. Therefore, we find out what the function 
is with only one application of the /-controlled-NOT gate. Classically, we would need to evaluate the 
function n times to find y. 

3 Learning which variables a general Boolean function depends on 

If f(x) is a general Boolean function, then when we measure \\\f out ) in the computational basis, we will 
obtain the label of one of the basis vectors v y , with which v/ has a nonzero overlap. The key to using this 
to solve the problem stated in the Introduction is the following fact: if f(x\,X2,- --x n ) is independent of 
the variable Xj, and y G {0, 1}" has the property that yj = \, then (v y |v/) = 0. In order to prove this we 
start by noting 



(>•>>•/> = ^ E (-D /M+ "' 



xe{0,l} n 



4^i---i(-i) /w+ *- y - (id 

v^v x x o 



Now look at the xj sum, 



£(_l)/W+*y = (_i)/W Y[ (-\) Xky * £(-1)^=0. (12) 

xi=0 k=l,tyj Xj=0 

This proves our result. What it implies is that if we use the Bernstein- Vazirani circuit with a Boolean 
function that is a junta and find an output vector \y) that has ones in a number of places, then the func- 
tion does depend on the variables corresponding to those places. If the function does not depend on a 
particular input variable, then the ra-qubit state \ \\f out ) will always have a in that position. 

It is important to note that the success probability to find the variables the function depends on is 
independent of the total number n of input variables. In general, the success probability for the quantum 
algorithm depends only on the form of the Boolean function that is being tested, that is, it depends on 
the number of significant variables, and the functional form of the Boolean function involving these 
significant variables. 
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3.1 Boolean functions depending on only two input variables 

In order to illustrate this, let us consider a simple example. Suppose that we know that our function is 
given by f(x\ ,X2, • • -x n ) — XjX/ci but we do not know j and k, i.e. we know that the Boolean function is 
the product of two of the variables, but we do not know which two. Our task is to find out which two. 
The vector |v/) corresponding to this function has a nonzero inner product with only four of the basis 
vectors |v y ). We must have y/ = for / / j,k, which leaves four possibilities, which we shall denote 
by \yoo), corresponding to yj = yt = 0, \yoi), corresponding to yj = and y^ = 1, etc. We find that the 
output of the Bernstein- Vazirani circuit in this case is 

Wont) = \{\ym) + boi) + bio) - |vn)). (13) 

If we measure in the computational basis, then we will obtain one of these basis vectors. If we obtain 
|yoo), we learn nothing, and the procedure has failed. This happens with a probability of 1/4. If we 
obtain either |joi) or bio), then we learn one of the variables, and if we obtain \yn), we obtain both. All 
of these outcomes have a probability of 1 /4, so that we learn at least one of the variables on which the 
function depends with a probability of 3/4. This probability is independent of how many input variables 
n there are in total. Classically, a possible procedure would be to initially set all of the variables equal 
to 1 , which would set the value of the function equal to 1 . One then changes the value of the variables, 
one at a time, to see which ones cause the value of the function to change. In order to learn on which 
variables the function depends, one would have to evaluate the function &{n) times. If n is large, the 
quantum procedure, though probabilistic, is more efficient. 

Let us now consider a somewhat more general example. We will still assume that our function only 
depends on two out of the n variables, xj and xu, but we will not assume the specific form of the function. 
We can express f(x\,X2, ■■■x n ) as 

f(x\,x 2 ,...x n ) =g(xj,x k ), (14) 

where g(xj,x^) is some Boolean function of two variables. Now, assuming that yi = for / / j,k, we 
have that 

i i 

( v f\ v y) = J L {-l) g{x ' A)+y ' x ' +ykXk . (15) 

xj,x k =0 

The right4iand side of the equation can only be 0, 1, or ±1 /2, and it will only be or 1 if f(x\ ,X2,-- -x n ) 
is one of the basis functions. Therefore, \y out ) is either one of the vectors \yi { i 2 ), or of the form 

Wout) = ^(±boo) ± boi) ± bio) ± bn))- (16) 

If f{x\ ,X2, ■ ■■Xn) is one of the basis functions, we succeed after one trial, however, we do not know this, 
and several trials in which we get the same answer will be necessary to confirm that we have one of 
the basis functions. In all other cases, we will fail, that is get no information about which variables the 
function depends on, with a probability of 1/4, so that after several trials we will, with high probability, 
know Xj and x^. 

3.2 Boolean functions depending on more than two input variables: an example 

Now let us see what happens if the function depends on more than two variables. We know that the 
quantum algorithm will always find the variables a function depends on, but that the success probability 
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for this will vary with the form of the Boolean function. Let us consider the case 



7/1 



f(x u x 2 ,...x n ) = Y[*j- (17) 

7=1 

The probability to identify which variables this function depends on would be the same also for other 
Boolean functions which are a product of any m out of the n variables. For vectors \v y ) such that yj = 
for j > m, we have 

(^M = it---I(-i) Mxi '- w) ) as) 

Z xi =0 x,„=0 

where 

m m 

h(xi,...x m ;y) =Y[ x j+Y* x j y i- ( 19 ^ 

7=1 7=1 

Now, if the product x\x 2 . ■ -x m were absent from the exponent in Eq. ( fT8l ), and if at least one of the yj ^ 0, 
then the sum would be zero. The product changes the sign of only one of the terms, so that we have 

< v /l*> = ±2^r C2Q) 

If yj = for j = 1,. . .n (we shall denote the vector corresponding to this y by |vo)), then without the 
product in the exponent all of the terms in the sum in Eq. (fT8l) would be 1. The presence of the product 
again changes only one term, so that 

(v/|vo> = (2D 

Note that since the failure probability is just |(v/|vo)| 2 , this implies that the failure probability grows 
with m. This is the "worst case scenario"; this type of Boolean function belongs to the class of functions 
for which the Bernstein-Vazirani algorithm has least probability to succeed in finding the variables it 
depends on, since a phase factor is added only to a single term. Nevertheless, the success probability is 
still independent of the total number n of input variables. 



4 Amplification of the success probability 

The desirable outcomes of the measurement of the output state \yfout) w& those with as many l's as 
possible, since a "1" in position i indicates that the Boolean function depends on input variable Xj. To 
further increase the success probability, it is possible to amplify components of \\\f out ) with a chosen 
number and above of l's. This procedure is based on Graver's quantum search algorithm. Graver's 
algorithm uses €?(N/M) queries for searching a database with iV elements, where M of these are solutions 
to the search problem J6l|7l. Classically, &(N/M) database queries are needed. 
Let us define the normalised states \a) and |j8) as 

L" 1 
v,|jc); A= (22) 

|/3) = b£v x \ X ); B = —^, (23) 
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where the prime ' indicates a sum over all x 6 {0, 1}" which contain k or more l's and " indicates a sum 
over the remaining x. The state | y ut) in terms of \a) and |j3) is 

11 

Wont) =j\ a ) + g\^) =cos-|a) + sin-|j8), (24) 



where cos = l/A = y £ v 2 . and sin = 1 /B = y £ vj. Repeated application of the operator 

G = //® n [/ / //®' l (2|0)(0| - l)H® n U f H® n O, (25) 
where the operator O produces phase factors —1 for components with k or more l's, gives 

G'\Yout) = cos l«> + sin (^"^) li 3 )- ( 26 ) 

after / applications. The optimal number of Grover iterations is given by the integer closest to 

, , arccos [sin (0/2)1 arccos,/y 

*w= — [ -^- LA = 2^= r < 27 > 

where y = £ v 2 . The leading term in the power series expansion of R(y) about y = is n/ (4-^/y). All 
higher order terms have a negative sign. Hence we have 

R< T7f- <28) 

and if y < 1, then 

R ~^rr <29) 

For this number of iterations, the final state contains the largest possible fraction of the component |j8). 
If the form of the Boolean function is known (e.g. that the Boolean function is of the form X{Xj, but not 
what i, j are), then it is possible to calculate y and the optimal number of Grover iterations for the chosen 
value of k. The smaller k is chosen, the larger y is, and the fewer Grover iterations are needed. If the 
form of the function is not known, then, just as for the usual Grover search algorithm, it is possible to 
estimate the optimal number of Grover steps [10]. This will require more queries of the function to be 
tested. It does, however, not necessarily mean that a significantly greater number of function queries is 
needed; this is the case for the example below. 



4.1 Amplification for a single term of order k 

As an example, let us consider the case where f(xi,Xz, - ■ -x n ) = Yl'j = \Xj, and suppose that we want to 
identify all variables this function depends on. As also pointed out before, the success probability would 
remain the same for any Boolean function which is a product of m input variables. From equation (20), 
we obtain y = 2~ 2m+2 , and consequently the optimal number of Grover iterations needed in order to 
obtain a high probability of identifying all input variables the function depends on is given by the integer 
closest to R = 71 2 m ~ 3 , which is 0(2 m ). Each iteration uses two queries of the Boolean function, so that 
the total number of function queries is roughly 2R = % 2 m ~ 2 , which is also &(2 m ). We point out that this 
number is independent of n, which is the total number of input variables. 
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If the Boolean function is a product of m of the input variables, but we do not know this, then we 
first need to estimate the optimal number of Grover iterations. It can be shown iPTOl that for a product of 
m input variables, the circuit for estimating the optimal number of Grover steps requires &(2 m ) function 
queries. In other words, if we are looking to amplify terms with m or more l's, that is, to find all 
variable the function depends on, then having to estimate the required number of Grover iterations does 
not change the order of how many function queries are needed in total. 

We can compare the success probability of the amplification strategy to the case where we run the 
unmodified Bernstein- Vazirani algorithm roughly 2R = % 2 m ~ 2 times (the number of runs is given by 
the integer closest to this number). In each round, the failure probability is (1 — 2~ m+1 ) 2 , so that the 
probability to fail in all rounds, learning none of the variables the function depends on, is approximately 
Pf = (1 — 2~ m+l ) n 2 . The probability to obtain at least one variable is therefore approximately l—p/, 
which approaches 1 — e~ n « 0.96 when m becomes large. On the other hand, the probability to never 
learn one particular variable x, that the function does depend on, in any of the 2R = n 2 m ~ 2 tries, is equal 
to 

(\ JZ2'"- 2 
I Kv/M 2 ] =(l-2- m+i r 2 "- 2 . (30) 
v,:y,=0 / 

This probability approaches e~ n l 2 rj 0.21 when m becomes large. For 2R function queries, there is 
therefore an appreciable probability for not learning at least one variable the function depends on when 
using the Bernstein- Vazirani algorithm without amplification. The amplified procedure is very likely to 
obtain all variables which the function depends on with a similar number of function queries. Amplitude 
amplification for terms with mi's has therefore improved the situation. 



5 Conclusions 

We have shown that the Bernstein- Vazirani algorithm may be used for testing which input variables an 
unknown Boolean function depends on. In a sense, this task is more general than distinguishing between 
linear Boolean functions, which is the task for which the Bernstein- Vazirani algorithm was originally 
devised. The success probability of finding variables a Boolean function depends on may be further 
enhanced by an amplification procedure based on Grover's search algorithm. 

The success probability for the presented quantum algorithm depends on the particular form of the 
Boolean function, but has the general property that it is independent of the total number of input variables. 
It shares this property with the algorithm presented in [9 ]. Nevertheless, a full comparison of the success 
probabilities of the different quantum and classical algorithms remains to be made. Other variations of 
the Bernstein- Vazirani algorithm may also be tailored for investigating Boolean functions of particular 
forms, and this will be the subject of further investigations. 
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